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Abstract 

For nonrelativistic Hamiltonians which are shape invariant, analytic expressions for the eigen- 
> ■ 

■ values and eigenvectors can be derived using the well known method of supersymmetric quantum 

mechanics. Most of these Hamiltonians also possess spectrum generating algebras and are hence 
solvable by an independent group theoretic method. In this paper, we demonstrate the equiva- 
qq ' lence of the two methods of solution by developing an algebraic framework for shape invariant 

Hamiltonians with a general change of parameters, which involves nonlinear extensions of Lie 
algebras. 



a: 

; 1 Introduction 

> ■ 

Supersymmetric quantum mechanics (SUSYQM) provides an elegant and useful prescription 
for obtaining closed analytic expressions both for the energy eigenvalues and cigcnfunctions of a 
large class of one dimensional problems. The main ingredients in SUSYQM are the supersymmetric 
partner Hamiltonians H- = A^A and H + = AA^ . The A and A'' operators used in this factorization 
are expressed in terms of the superpotential W as follows: 

A(x, a) = — + W(x, a) ; A\x,a) = — ^- + W(x,a). (1) 
ax ax 

Here, W is a real function of x and a is a parameter (or a set of parameters), which plays an 
important role in the approach of this paper. An interesting feature of SUSYQM is that for a shape 
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invariant system |2|, i.e. a system satisfying an integrability condition 

W\x, oo) + ^^9l = W\x, ai ) - d -^pl + R(a ) ; a, = /(a ) , (2) 

the entire spectrum can be determined algebraically without ever referring to underlying differential 
equations jy. 

Several of these exactly solvable systems are also known to possess what is generally referred to 
as a spectrum generating algebra (SGA) || M. The Hamiltonian of these systems can be written as 
a linear or quadratic function of an underlying algebra, and all the quantum states of these systems 
can be determined by group theoretic methods. 

One may naturally ask the question whether there is any connection between a general shape 
invariance condition and a spectrum generating algebra. In this paper we address the equivalence 
between the two approaches, considering a large class of change of parameters, including translations, 
scalings, projective transformations, as well as more complicated functions /(ao). 

In sec. 2, we start with a general shape invariant model. We make use of the operators A 
and A' to construct a three generator algebra. In particular, the shape invariance condition plays 
a crucial role in closing the algebra, which turns out to be either so(2, 1) or a deformation of it. 
In sec. 3, several examples are presented. In particular, we discuss shape invariant potentials 
generated by a change of parameters corresponding to translation a± = ag + k and pure scaling 
a± = qcio,q = constant (0 < q < 1). For the case of scaling, we find that the associated potential 
algebra is a nonlinear deformation of su(2). We also describe potential algebraic structure of cyclic 
potentials |^| described as a series of shape invariant potentials which repeats after a cycle of k 
iterations. And finally, we discuss the potential algebra of Natanzon potentials and show that all 
translational shape invariant potentials can be generated from them. 



2 The Algebraic Shape Invariant Model 

To begin the construction of the operator algebra, let us express the shape invariance condition eq. 
(H) in terms of A and : 

A(x, ao)A\x, ao) — A^(x, ai)A(x, a±) = R(ao) ■ (3) 

This relation resembles a commutator structure. To obtain a closed su(2)-like algebra, we introduce 
an auxiliary variable (f> and define the following operators 

J+ = e i **Ai(x, X (id4,)) , J-=A(x,x(idt))e- i ^ , (4) 

where p is an arbitrary real constant and x is an arbitrary, real function. The operators A{x, x(id(/>)) 
and A^(x, x{i®(f>)) are obtained from eq. ([!]) with the substitution a — » x(*^)- This general- 
ization is analogous to the familiar spherical coordinate separation of variables scheme, in which 
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9^$(<^>) — > constant x $(<^>); in this case the constant eigenvalue is ao- From eq. (Q), one obtains 

[J+,J_] =e i ^t( s>x (0 ))^( XjX (^)) e -iP*_ i4 ( a;)X ( i ^)) i4 t( a;)X ( i ^)) . ( 5 ) 

Eq. (|5|) can be easily cast into the following form 

[J+, J_] = - {A(ar, x (i^))At( a; , x( ^)) _ A t (Xj x(i ^ +p))A( a; , x (i^ + p ))} . (6) 

At this point if we judiciously choose a function x{^4>) such that x(^4> +p) = /[x(*^)L the r.h.s. 
of eq. (^) can be simplified using shape invariance condition 

Afox^JAtfoxtify)) _ At(x,x(^ +p))A(i )X (<^ + P)) = RUiid*)) , (7) 

where we have identified 

<*o x(ify) ; «i = /(cto) /[x(ify)] = . (8) 

The last step in our construction is to define the operator J3 as J3 = — . As a consequence, we 
obtain a deformed Lie algebra whose generators J+, J_ and J3 satisfy the commutation relations 

[J 3 ,J±] = ± J± ; [J+,J-]=e(J 3 ) , (9) 

where £(Js) = — R(x(id<f>)) defines the deformation. Thus we see that shape invariance condition 
plays an indispensible role in the closing of this algebra. 

Depending on the choice of the x function in eq. (ih, we have different reparametrizations 
corresponding to several models. For example we have 

1. translational models: a\ — ao + p for x( z ) = z ( m these models if R is a linear function 
of J3 the algebra turns out to be so(2, 1) or so(3) [[?] ; a similar conclusion was reached by 
BalantekinQ by using a somewhat different method;) 

2. scaling models: a\ = e p a,Q = qao for x( z ) = e z , 



3. cyclic models: a, = ^±§ , for = ^£)^H^)M^ , 

where Ai ; 2 are solutions of the equation (x — a)(x — 8) — (3"f = and B(z) is an arbitrary periodic 
function of z with period p. 

Other changes of parameters follow from more complicated choices for x( z )- For example, if one 
takes x( z ) = e e j one §> e ^ s the change of parameters a\ = Oq. 

Note that the quantity J+ J_ corresponds to the Hamiltonian 

H-(x,id$+p) = A^(x,x(i9<t,+p))A(x,x(id^,+p)) . (10) 
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To find the energy spectrum of the Hamiltonian H- of cq. fllOj), we first construct the unitary 
representations of the deformed Lie algebra defined by eqs. (^). The technique proceeds as follows 
||. Define, up to an additive constant, a function g{Js) such that 

£(J S ) = g(J 3 ) - g(J 3 - 1) ■ (11) 

The Casimir of this algebra is then given by C2 = J~ J+ +g{ ^3)- It is known that in a basis in which 
J3 and C2 are diagonal, J + and J_ play the role of raising and lowering operators, respectively. 
Operating on an arbitrary state \h) we have 

J a \h) = h\h) , 

J_|ft) =a(h) \h-l) , 

J+\h) =a*(h + l) \h + l) . (12) 

Using eqs. ([}]) and (|l^) we obtain 

\a(h)\ 2 -\a(h + l)\ 2 =g(h)-g(h-l) . (13) 

The profile of g(h) determines the dimension of the unitary representation. For example, let us 
consider the two cases presented in fig. One obtains finite dimensional representations fig. la, 




(a) (b) 

Figure 1: Generic behaviors of g(h). 

by starting from a point on the g(h) vs. h graph corresponding to h — h m i n , and moving in integer 
steps parallel to the /i-axis till the point corresponding to h — h max . At the end points, a(h m i n ) = 
a(h ma x + 1) = 0, and we get a finite representation. (This is the case of su(2) for example, where 
g(h) is given by the parabola h(h + 1).) If g(h) is decreasing monotonically, fig. lb, there exists 
only one end point at h — h m i n . Starting from h m in the value of h can be increased in integer 
steps till infinity. In this case we have an infinite dimensional representation. As in the finite case, 
h m in labels the representation. The difference is that here h m i n takes continuous values. Similar 
arguments apply for a monotonically increasing function g(h). 
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Having the representation of the algebra associated with a characteristic model, we obtain (using 
eq. ( |Io|Jl^ ) ) the complete spectrum of the system. To illustrate how this mechanism works, we 
investigate few examples in the next section. 

3 Examples 

3.1 Self-Similar Potentials 

The first example is for a scaling change of parameters a\ — qa Q . Consider the simple choice 
R(a>o) — rido, where r\ is a constant. This choice generates self-similar potentials studied in refs. 
0,0. In this case, eqs. (KJ) become: 

[J 3) J±] = ± J± ; [J+, J_] = £(J 3 ) = -n exp(- P J 3 ) , (14) 

which is a deformation of the standard so(2, 1) Lie algebra. 
For this case, from eqs. (Oh and (O) one gets 



9(h) = ^p-re-v h = --^q- h ; q = <* . (15) 
e p — 1 1 — q 

Note that for scaling problems jy}, one requires < q < 1, which leads to p < 0. From the 
monotonically decreasing profile of the function g(h), it follows that the unitary representations of 
this algebra are infinite dimensional. If we label the lowest weight state of the operator J3 by h m i n , 
then a(h m i n ) = 0. Without loss of generality we can choose the coefficients a(h) to be real. Then 
one obtains from ( |l3| ) for an arbitrary h = h m in + ?l i n = 0, 1, 2, . . . 

a 2 (h) =g(h-n-l)- g(h - 1) = n ^ q 1 ^ 1 . (16) 

9-1 

The spectrum of the Hamiltonian H_(x, a±) is given by 



H-\h) = a 2 (h)\h) = n ^— - q 1 -' 1 ^) . (17) 



Therefore, the eigenenergies are 



a n — 1 

E n (h) = naih)^—^ ; a(h) = q 1 ^ . (18) 

To compare the above spectrum obtained using a group theoretic method with the results obtained 
from SUSYQM Jll[], we go to the x-representation. Here \h) oc e. lph ^iph mirl ,n(x) and hence, the 
Schrodinger equation for the Hamiltonian if_ reads 

|-^2 + W2 ( x > e ^ +P ) - e W " +p ) - E J e^ /l ^ ml „,r l (x) = , 

or 

~ + W 2 (x, a(h)) - W'(x, a(h)) - e\ ^ mm ,„(i) = , (19) 
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which is exactly the Schrodinger equation appearing in ref. |pl| , with eigenenergies given by eq. ([l8|). 
The elegant correspondence that exists between potential algebra and supersymmetric quantum 
mechanics for shape invariant potentials is further described in ref. ]T^ | . 
For a more general case, we assume i?(ao) = Sj=i ^' a o- I n this case 

oo „ 

7=1 

and therefore one gets 



a 2 (^) = <?(/>- n-1) 

1 - q> n 



where ctj(h) = i^e -3 '^ -1 ). These results agree with those obtained in ref. Ell . 
3.2 Cyclic Potentials 

Let us consider a particular change of parameters given by the following cycle (or chain) : 

ao, cti = f(a ), a 2 = /(ai) , . . ., a fc _i = /(a fe _ 2 ), a fc = f(a k -i) = a , (22) 

and choose B(a,i) = a, = Wj. This choice generates cyclic potentials studied in ref. j^j. 

Cyclic potentials form a series of shape invariant potentials; the series repeats after a cycle of k 
iterations. In fig. |^ we show the first potential V(x, ao) from a 3-chain (k = 3) of cyclic potentials, 
corresponding to uoq — 0.15, ui% — 0.25,^2 = 0.60. 

Such potentials have an infinite number of periodically spaced eigenvalues. More precisely, the 
level spacings are given by w 0) wi, ■ ■ ■ ,<Jk-i,wa, coi, ■ ■ • , u)k-i,Wo, 

In order to generate the change of parameters ( p2| ) the function / should satisfy /(/(•• -f(x). . .)) = 
f h (x) = x. The projective map 

M = ^| , (23) 
IV + o 

with specific constraints on the parameters a,/3,7, S, satisfies such a condition ||. 

The next step is to identify the Lie algebra behind this model. For this, we need to find the 
function x satisfying the equation 

X(* + P) = /(x(z))= a *Ktf ■ ( 24 ) 
It is a difference equation and its general solution is given by 

x[z) = ^-Srf' + ^-Srf'Bjz) 



7 X^ + X; /p B(z) 
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Figure 2: First potential V(x, ao) from a 3-chain (k = 3). 

where Ai_2 are solutions of the equation (x — a)(x — 5) — (3"f = . For simplicity B(z) can be chosen 
to be an arbitrary constant. Plugging this expression in eqs. (j^) we obtain: 



[J+, J-] = «*) = _l^x-^ + B(A a -^- 



(26) 



A\ t- x^ 2 

Applying our standard procedure to find the spectrum of the Hamiltonian H- = J+J- we find 
that the ground state is at zero energy; the next [k — 1) eigenvalues are Ei = Ylj=o > I ~ 
0, 1, . . . , (k — 2) , and all other eigenvalues arc obtained by adding arbitrary multiples of the quantity 
Qft = ujq + lui + ■ ■ ■ + 0Jk-i- This result is in complete agreement with H. 



3.3 Natanzon Potentials 

In sec. 2, we noted that for SIP's with translationally related parameters (i.e. a± = ao + 1), the 
shape invariance condition helps in closing the algebra to the familiar so(3) or so(2, 1), provided 
the R(ao) was linear in ao Q. Several SIP's belong to this category; among them are the Morse, 
Scarf I, Scarf II, and generalized Poschl- Teller potentials. However, there are many important 
SIP's (e.g., Coulomb), whose associated R(ao)'s are not linear in ao- Our method of the previous 
section would lead to deformed potential algebras for these systems. While we now know how to 
get representations of such algebras, in this section we shall take a different approach. We choose 
to generalize the structure of operators J± such that their algebra still remains linear. In fact, in 
this section, we generate shape invariant potentials from an underlying potential algebra instead of 
showing algebraic structure hidden in a shape invariant system. 
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Alhassid ct al. || had shown that the algebra associated with the general potential of the Natanzon 
class is so(2,2). The Schrodinger equation for these potentials reduce in general to the hyperge- 
ometric equation. We show below that a further constraint generates all SIP's with translational 
change of parameters. For the sake of completeness we will briefly examine the properties of so(2, 2) 
algebra in this section, and show its connection to the Natanzon potentials ||. We then conjecture 
an additional constraint that would render them shape invariant. We find that this conjecture indeed 
produces all known SIP of the translational type. We shall find in fact that the subset of Natanzon 
potentials associated with the translational (additive) SIP's has the simpler so(2, 1) algebra. 

We begin by describing Alhassid et al.'s realization of the so(2, 2) algebra in terms of differential 
operators. For consistency, we use the formalism and the notations of refs. ||. 

The differential realization can be written explicitly as 



I 

2 6 



±i(<P+e 



d 



Tq- + tanhx M<9</,) + cothx (-«<%) 



-- (dt + d e ) ; 



1 



>±i(0-0) 



d 

+ tanhx {-id<j>) + cothx (+id 9 ) 



A± = Ai± A 2 
A 3 

B± = B\ ± B 2 
B 3 

The so(2, 1) algebra obeyed by these operators is 

[A 3 ,A±} = ±A ± , [A + ,A_] = -2A 3 , 
and a similar one for the £Ts. The Casimir operator Ci is given by 



(27) 



-^{d^-dg) 



Co 



2 (A 



A+A^ - A 3 



2 (Bl 



B+B- — B 3 



& 



d 



-—r + (tanhx + cothx) t; 1" scc h X (~ id<i>) ~ cosech x ( — ids 



(28) 



_ t>x ox 
Operators ^3, B 3 and C2 can be simultaneously diagonalized, and their actions on their common 
eigenstate are given by 



C 2 |w,mi,TO 2 ) = uj(u + 2) \u,m 1 ,m 2 ) 
A 3 \u!,m 1 ,m 2 ) = mi \u>,mi,m 2 ) ; 
B 3 \Lo,m ll m 2 ) = m 2 \u},m 1 ,m 2 ) ■ 



(29) 



(It is worth mentioning at this point that the Casimir operator given above is indeed self-adjoint 
with respect to a measure sinhxcoshx^X^^-) 

Now we shall briefly describe a general Natanzon potential and show its connection to the above 
Casimir operator. A general Natanzon potential U (r) is implicitly defined by p] 



U[z(r)] = -f^-z)+^(l- Z ) + h lZ _ 1 { ^ r} 



(30) 
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with Q(z) quadratic in z: Q(z) = az 2 + boz + co = a(l — z) 2 — — z) + a and /, ho, hi, a, bo, b\, Co, c\ 
are constants. The Schwarzian derivative {z, r} is defined by 

(31) 



, , d 3 z/dr 3 
{z,r} = 



dz/dr 



d 2 z/dr 2 



dz/dr 



The relationship between variables z (0 < z < 1) and r is implicitly given by 

dz\ _ 2z(l - z) 



(32) 



To connect the Casimir operator Ci of the so(2, 2) algebra [eq. (|28h] to the general Natanzon 
potential, we perform a similarity transformation on Ci by a function F and then follow that up by 
an appropriate change of variable x = <?(?*)• It has been shown JL3| that to turn Ci into the form of 



a Schrodinger Hamiltonian, one needs to choose F 
EQ 



(sinh( 



h(2g) V 



and choose z = tanh g. Then 



U(z(r)) = 



-1 + ^-z 



-z 2 \ 

4^ \ 



z(l - z) {-id^f + (1 - z) {-idgf 1 

-Q 2 {z ' r} 



We have used 



- (aE - ~ + (-^) 2 J z(l - z) + - J + (-z9 e ) 2 j ( I - : ) 

+ ((a + b + co)E - 1)] /Q(z) - \{z,r) . 



, dg dg dz 



dr dz dr 2y / z(l — z) Q 



2z(l - z) _ IT 



g = tanh 1 \fz and from eq. (|32]). Now, with the following identification 

/ = aE - 7 - + H9 ) 2 , 

7 2 
h = coE -- + (-ide) , 

hi = (a + b + c )E-l , 



(33) 



(34) 



the potential of eq. ( |33| ) indeed has the form of a general Natanzon potential [eq. (J30|) ] - Further 
details are shown in ref. Finally, we are ready to explicitly demonstrate the connection 

between the Natanzon potential algebra and shape invariant potentials of supersymmetric quantum 
mechanics. 

We now note that the similarity transformation can be rewritten: using g = tanh -1 ^fz, g' = 
and eq. (|32"1), we find ( EHMHgl J — j*.. At this point we go back to the operators A± [eq. (|27j)1 and 

/4. This 

\ » / * Z 

transformation carries operators A± to 



rinh(2g) \ 

9' y 

1 

ask how they transform under the similarity transformation given by F ~ ^ sin M 2 fl) ^ 



A, 



<i 1 dz' 1 dz 
c?X 2z' dx 2z 



+ tanhx {—id(j>) + cothx {-ide) 



(35) 
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Except for the expression {^j^§^ — 573^)1 this looks very much like eq. (p7f), which arc in fact 
A± of the shape invariant Poschl- Teller potential 0. Thus, if f^F^r — S3x) were ^° ^ e a l mear 
combination of tanh% and coth%, operators A± could be cast in a form similar to the operators A± 
of eq. (|27]), and we would get A±'b that generate shape invariant Hamiltonians. 
Hence to get shape invariant potentials, we require, 



1 dz' 1 dz 
2z'~d~X ~ 2z~d~x 



atanhx + /3cothx- (36) 



This leads to 

z' = z 1+(i (l - z)- a - p , (37) 
which is the second constraint on the relationship between variables z and r. Since these variables 



are already constrained by eq. (32), only a handful of solutions would be compatible with both 
restrictions. The z(r)'s that are compatible with both eqs. ( |32"| ) and (|37|)are given by 

_ z) -«-/» = 2£g^*) | (38 ) 

vQ( z ) 

where Q(z) is a quadratic function of z. After some computation, we find that there is only a finite 
number of values of a, f3 which satisfy eq. (p8|). These values are listed in Table 1, and they exhaust 
all known shape invariant potentials that lead to the hypergeometric equation. 

Furthermore, while the potential algebra of a general Natanzon system is so(2, 2), and requires 
two sets of raising and lowering operators A± and B±, all translational shape invariant potentials 
need only one such set. For all SIPs of Table 4.1 of ref. |Q, one finds that all partner potentials 
are connected by change of just one independent parameter (although other parameters which don't 
change are also present.) Thus there is a series of potentials that only differ in one parameter. From 
the potential algebra perspective, all these potentials differ only by the eigenvalue of an operator 
that is a linear combination of A3 and B3, and all are characterized by a common eigenvalue of C2. 
Thus, these shape invariant potentials can be associated with a so(2, 1) potential algebra generated 
by operators A+, A^ and the same linear combination of A3 and B3. 

Conclusion: In this paper, we have explored the reasons underlying the integrability of shape 
invariant Hamiltonians and shown that such systems naturally admit an algebraic structure known as 
potential algebra. We have derived these algebras for shape invariant systems with translational and 
scaling type change of parameters, as well as for cyclic potentials. In general, one finds deformations 
of the so(2, 1) Lie algebra. Our approach links the group theoretic and supersymmetric quantum 
mechanics approaches for treating shape invariant potentials. 

A.G. acknowledges a research leave and a grant from Loyola University Chicago which made his 
involvement in this work possible. Partial financial support from the U.S. Department of Energy is 
gratefully acknowledged. 
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a 


(3 


z(r) 


Superpotential 


Potential 








z = e~ r 


m! coth | + m 2 


Eckart 





i 

2 


• 2 r 

z = sin 2 


micosec r + to 2 cot r 


Gen. Poschl- Teller 
trigonometric 





-1 


z = 1 - e- r 


mi coth | + m 2 


Eckart 


l 

2 





z = sech 2 | 


TOiCosechr + to 2 cothr 


Poschl-Teller II 


1 

2 


1 

2 


z = tanh 2 | 


mi tanh § + m 2 coth | 


Gen. Poschl-Teller 


-1 





z = 1 + tanh | 


mi tanh 1+7712 


Rosen Morse 



Table 1: Shows all allowed value of a, (3 and the superpotentials that they generate. 



12 



